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The heat transport of the spin-gapped material (CH 3 ) 2 CHNH 3 CuCl 3 (IPA-CuCls), a candidate 
quantum magnet with Bose-Einstein condensation (BEC), is studied at ultra-low temperatures and 
in high magnetic fields. Due to the presence of the spin gap, the zero-field thermal conductivity (k) 
is purely phononic and shows a ballistic behavior at T < IK. When the gap is closed by magnetic 
field a,t H = Hd, where a long-range antiferromanetic (AF) order of Cu^+ moments is developed, 
the magnons contribute significantly to heat transport and exhibit a ballistic T^ behavior at T < 

600 mK. In addition, the low-T k{H) isotherms show sharp peaks at Hd, which indicates a gap 
re-opening in the AF state {H > Hd) and demonstrates limited applicability of the BEC model to 
IPA-CuCls. 

PACS numbers: 66.70.-f, 75.47.-m, 75.50.-y 

I. INTRODUCTION 


Low-temperature thermal conductivity (k) can probe 
the transport properties of various elementary excitations 
in solidsAi^ At very low temperatures, the phonon trans¬ 
port is known to exhibit a ballistic behavior since all the 
microscopic scatterings are smeared out and phonons are 
only scattered by the sample surface or boundary. This is 
the so-called boundary scattering limit and the n shows 
a simple dependence^— The magnon excitations of 
an antiferromagnetically ordered material are known to 
have the same statistic law as the phonons. In addition, 
since the low-energy antiferromagnetic (AF) magnons 
also have linear dispersion, they are able to show the 
same dependence of thermal conductivity at very low 
temperatures—"— However, this ballistic transport of AF 
magnons has been rarely observed— One reason is that 
the low-energy magnons and phonons can easily couple to 
each other—""— More serious is that the acoustic magnons 
are almost always gapped in the real antiferromagnets. 
In this regard, it is possible to switch on the magnon 
heat transport by applying magnetic field to close the 
anisotropy gap— For example, in the AF insulator 
Nd 2 Cu 04 , a 10 T-field-induced increase of k was found to 
show a behavior at ultra-low temperatures, which was 
discussed to be the first observation of the AF magnon 
ballistic transport/ However, this explanation is not well 
grounded since the anisotropy gap can be closed only at 


the spin re-orientation field (much lower than 10 T for 
Nd 2 Cu 04 ) and it will be re-opened at higher field— 

It might be easier to probe the ballistic transport of 
the AF magnons in those quantum magnets that can ex¬ 
hibit the Bose-Einstein condensation (BEC)— BEC 
denotes a collective occupation of bosons to the lowest 
single-particle state when temperature approaches zero. 
In some spin-gapped quantum magnets, the XY-type AF 
state induced by the magnetic field that closes the spin 
gap can be described as a BEC state. The spin Hamilto¬ 
nian of these quantum magnets must contain the U{1) 
symmetry, which requires a continuous uniaxial sym¬ 
metry. In field-induced BEC state, the U{1) symme¬ 
try spontaneously gets broken and, as a consequence, 
a gapless Goldstone mode is acquired. If this kind of 
low-energy AF magnons are not strongly coupled with 
phonons, they are likely to contribute to transporting 
heat. 

Probing the magnon heat transport of BEC materi¬ 
als have been tried for the organic compound NiCl 2 - 
4SC(NH2)2 (DTN), an S' = 1 chain system /°i^^ and 
the oxide Ba3Mn208— ^ It has been found that along 
the chain direction of DTN, the k shows a distinct en¬ 
hancement when increasing field across the BEC phase 
boundary— However, the reported data could not sep¬ 
arate the magnon thermal conductivity from the total 
K, mainly because the measurements have not been car¬ 
ried out at low enough temperatures— In the case of 
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FIG. 1: (Color online) Schematic H — T phase dia¬ 
gram of IPA-CuCls according to the magnetization, spe¬ 
cific heat, magnetocaloric effect, and neutron scattering 
measurements.— PM represents the low-field quantum 
paramagnetic state and AFM represents the field-induced an- 
tiferromagnetically ordered state, which has been discussed to 
be a possible BEG state. The solid line is the phase boundary 
(Hci) of the above two states. Note that the critical fields 
10 T at T —>■ 0) are slightly different for different field direc¬ 
tions and from different experiments. Inset: Schematic plot 
of the ac-plane structure, including the magnetic Cu^’*' (big 
and red circles) and the bridging Cl“ (small and cyan circles) 
ions. Other atoms are omitted for clarity. 


Ba3Mn208, the magnons in the BEG state are scattering 
phonons rather than transporting heati^ In this work, 
we choose another BEG candidate (GH3)2GHNH3GuGl3 
(IPA-GuGl3)2^-.2a to study the thermal conductivity. 
IPA-GUGI3 crystallizes in a triclinic structure and the 
Gu^’*' spins {S = 1/2) form ladders along the a axis, 
with rungs along the c axis, as shown in Fig. 1. The 
zero-field ground state is quantum paramagnetic with a 
spin gap of 1.17 meV. When the magnetic field closes the 
gap at fioHci ~ 10 T, an AF state is developed. This has 
been proposed as a BEG state since the neutron scatter¬ 
ing indicated a gapless mode— >2^ Here, we study the heat 
transport of IPA-GuGls single crystal at very low temper¬ 
atures down to several tens of millikelvin. At H = Hd, 
a magnon thermal conductivity is observed at T < 
600 mK, which is the clearest experimental evidence of 
the ballistic magnon heat transport in the AF state till 
now. However, the low-T k{H) isotherms indicate that 
the spin-gap is closed only at Hd- The re-opening of a 
small gap at H > Hd demonstrates that the BEG model 
is not strictly applicable to IPA-GuGls. 


II. EXPERIMENTS 

IPA-GuGls single crystals were grown by a slow evap¬ 
oration of ethanol solution mixed with GuGl2, iso¬ 
propyl amine, and concentrated HGl in an appropriate 
proportion— The as-grown crystals are dark brown with 
size uptol5x5x0.6 mm^, with the length and width 
nearly along the c and a axis, respectively. The specific 
heat of an IPA-GUGI3 single crystal was measured by 
the relaxation method in the temperature range from 2 
to 30 K using a commercial physical property measure¬ 
ment system (PPMS, Quantum Design). For anisotropic 
K measurements, the long-bar shaped samples were cut 
from the as-grown crystals along either the a or c axis, re¬ 
spectively. The K was measured using a “one heater, two 
thermometers” technique in a ^He refrigerator at 300 mK 
<T < 30 K and a ^He-^He dilution refrigerator at 70 mK 
< r < 1 K, equipped with a 14 T magneti^ '^^i^^ In these 
measurements, the magnetic field is parallel to the heat 
current (Jr), which is along the lengths of the samples. 
It should be noted that the crystals are so fragile that 
cutting and polishing can easily produce some damages 
inside the samples; furthermore, some samples were even 
broken after cooling and warming cycle. Thus, the data 
presented in this paper were taken on several different 
samples. 


III. RESULTS 


Figure 2(a) shows the specific heat of an IPA-GUGI3 
single crystal with 2 < T < 30 K. Our data are con¬ 
sistent with the reported data— It is notable that the 
temperature dependence strongly deviates from the well- 
known T^-behavior of phonon specific heat. A fitting of 
low-T data to Cp = is shown by a dotted line in Fig. 
2(a). Furthermore, the data cannot be well fitted by a 
more complicated formula of phononic specific heat, 

Cp = I3T^ + + PrT\ ( 1 ) 


which is the low-frequency expansion of the Debye func¬ 
tion with /3, and P^ the T-independent coefficients— 
The reason is mainly due to the contribution of mag¬ 
netic excitations. As suggested by Ref. [H, since the 
existence of a spin gap of magnetic excitations, the low- 
T specific heat can be described as a T^ phononic term 
plus a Schottky term, that is. 


c = pT^ + in 


A e.^lkBT 

l^J (l-b(l/3)e‘^/'=Br)2’ 


( 2 ) 


where A is the spin gap separating the singlet ground 
state and the triplet first-excitation state, A is an ad¬ 
justing parameter. The factor 1/3 is related to the triple 
degeneracy of the first-excitation state— This formula 
can fit the low-T specific heat data well with the param¬ 
eters P = 4.81 X 10“^ J/K"'^mol, A/ks = 22.1 K and 
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FIG. 2: (Color online) (a) Low-temperature specific heat of 
an IPA-CuCla single crystal in zero field. The solid line is a 
fitting to the data using Eq. © with adjustable parameters 
P = 4.81 X 10"® J/K'^mol, A/Icb = 22.1 K and A = 0.88 
J/Kmol. The dashed line is a different fitting using Eq. @ 
with the fixed parameter A/fcs = 13.6 K and adjustable pa¬ 
rameters P = 7.18 X 10“^ J/K^mol and A = 0.19 J/Kmol. 
The dotted line is a simple T® fitting with coefficient of 7.92 
X 10“^ J/K'^mol. The inset shows data in a broader temper¬ 
ature range from 2 to 30 K. (b) Temperature dependence of 
thermal conductivity of an IPA-CuCla single crystal in zero 
field and at 300 mK - 30 K. The heat current is applied along 
the c axis. The dimension of this sample is 4.6 x 1.74 x 0.88 
mm^. The dashed line indicates T® temperature dependence. 
Inset: the low-temperature data in a linear plot for k vs T®. 
The thin line is a fitting to «: = bT^ with the parameter b = 
0.32 W/K*m for T < 700 mK. 


A = 0.88 J/Kmol, as shown in Fig. 2(a). These param¬ 
eters are also consistent with the results in the earlier 
literature;^ However, the fitting parameter A is rather 
different from the size of the spin gap (13.6 K) obtained 
from neutron measurements;^ If we fit the data using 
formula (©) with a fixed parameter A/ks = 13.6 K, the 
fitting is much worse, as shown by the dashed line in Fig. 
2(a). This discrepancy should be due to the complexity 
of low-T specific heat. Some other factors, such as the 
spin disorders or magnetic/nonmagnetic impurities, may 
give small additional contributions to the specific heat. 


Figure 2(b) shows the temperature dependence of k of 
an IPA-CuCla single crystal in zero field with the heat 
current along the c axis. Due to the spin gap of 1.17 meV 
(13.6 K) in the ground state of IPA-CuCla^i^i the mag¬ 
netic excitations can hardly be thermally excited at very 
low temperatures. Therefore, the low-T k{T) is the pure 
phonon conductivity. Indeed, the low-T k(T) behaves 
as a common insulator, with a phonon peak locating at 
about 6 K. The magnitude of peak exceeds 10 W/Km and 
is comparable with some other organic low-dimensional 
magnets 

A notable phenomenon is that the k(T) follows a per¬ 
fect T^ dependence at T < 700 mK. In principle, the 
phonon thermal conductivity can be expressed as a ki¬ 
netic formula Kp = ^Cvplp, in which C = PT^ is the low- 
T specific heat, Vp is the averaged sound velocity and is 
nearly T-independent at low temperatures, and Ip is the 
mean free path of phonons.— With decreasing tempera¬ 
ture, the microscopic scattering of phonons are gradu¬ 
ally smeared out and the Ip increases continuously untill 
reaching the averaged sample width W = 2\J where 
A is the cross-section area of sample;^ This boundary 
scattering limit of phonons can be achieved only at very 
low temperatures and the T-dependence of Kp is the same 
as the T^ law of the specific heat;i>^ The well-known T^ 
ballistic behavior of phonons, however, has been rarely 
observed in the transition-metal compounds, including 
the high-Tc cuprates, the multiferroic manganites and 
the low-dimensional quantum magnets.— 

In the boundary scattering limit, the phonon thermal 
conductivity of an isotropic system is given by^ 


27r2 Ip 



( 3 ) 


The averaged phonon velocity Vp can be extracted from 
the specific-heat coefficient using the relations /3 = 
and Sd = where Od is the 

Deybe temperature, N is the number of molecules per 
mole and each molecule comprises s atoms, V is the 
volume of crystal and R is the universal gas constant. 
For phonons in a three-dimensional (3D) lattice, the 
anisotropy of phonon velocity is usually not very strong 
and formula (©) can describe the Kp rather well. 

The present low-T data is fitted well to k = bT^ with 
b = 0.32 W/K^m, as shown in the inset to Fig. 2(b). 
Thus, with the parameter /3 = 4.81 x 10“^ J/K'^mol, the 
mean free path is calculated to be 0.017 mm. This value 
is much smaller than the geometry size of this sample, 
W = 1.40 mm. Since the T^ behavior is a robust signa¬ 
ture of the boundary scattering limit, it seems that the 
sample has some small cracks that act as the boundaries 
of phonon transport. As mentioned above, it is likely that 
this kind of cracks are not intrinsic but are produced in 
the cutting, polishing and cooling processes. 

Figure 3 shows the magnetic-field dependencies of Ka 
and Kc at low temperatures for two IPA-CuCla single 
crystals. The data exhibit two remarkable features. 
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FIG. 3: (Color online) Magnetic-field dependencies of the Ka 
(a) and «c (b) of two IPA-CuCls crystals at subkelvin tem¬ 
peratures. Since the magnetic fields are applied along the 
direction of the heat current, the demagnetization effect is 
negligibly small for these long-bar shaped samples. 


First, the k are field independent at low fields but ex¬ 
hibit a sharp peak at ^oHci ^ 9.95 and 9.75 T for H \\ a 
and II c, respectively, particularly at very low tem¬ 
peratures. It is noted that at T —>• 0 the peak positions 
are coincided with the transition field of the field-induced 
AF stateSince the magnons can be easily ex¬ 
cited when the spin gap is closed at Hd, the increase 
of the K at Hd is a direct evidence of the magnon heat 
transport. Note that the small anisotropy (9.95/9.75 = 
1 .02) of critical fields between H \\ a and || c is due to 
the anisotropy of g values. The electron spin resonance 
(ESR) measurements showed that the g values along 
the a and c axes are 2.05 and 2.11, respectively;^ This 
anisotropy (2.11/2.05 = 1.03) explains the anisotropy of 
Hd- Second, the peak feature demonstrates that the 
gap is closed only at Hd and is re-opened above Hd, 
which results in the vanishing of magnon transport in 
high fields. 

To probe the T-dependence of the magnon heat trans¬ 
port, the Ka(T) and Kc(T) at Hd of these two samples 
are measured and compared with their zero-field data, 
as shown in Fig. 4. Both the Ka and Kc in H = 0 fol¬ 




FIG. 4: (Color online) (a) Temperature dependencies of the 
Ka in zero field and 9.95 T (|| a), (b) Temperature depen¬ 
dencies of the Kc in zero field and 9.75 T (|| c). The data are 
taken on two IPA-CuCls crystals, of which the k{H) data are 
shown in Fig. 3. The dimensions of the Ka and Kc samples are 
4.5 X 1.29 X 0.55 mm® and 3.91 x 1.54 x 0.55 mm®, respectively. 
(c,d) Data shown with k vs T® plot. The solid lines are linear 
fittings to the experimental data. The dot dashed lines denote 
the magnon thermal conductivity. Km,a = 2.08T® (W/Km) 
and Km,c = 0.29r® (W/Km), obtained by subtracting the 
zero-field data from the k{T) curves in the critical fields. 
The dashed lines show the calculated curves Km,a = 2.27r® 
(W/Km) and Km,c = 0.26r® (W/Km), using formula ((5| with 
appropriate parameters, that are closest to the experimental 
data of Km,a and Km,c- 


low the K = bT^ dependence at low temperatures, with 
b = 0.82 W/K"^m for Ka and 0.50 W/K'^m for Kc, re¬ 
spectively. The phonon mean free paths in the ballistic 
regime are calculated to be 0.043 and 0.026 mm for the 
Ka and Kc samples, respectively, which are also smaller 
than the sample sizes {W = 0.94 and 1.03 mm). All these 
are consistent with the result of another sample shown 
in Fig. 2. 

The most important finding is that the Ka {kc) at 9.95 
(9.75) T also follow the T® dependence at T < 600 (700) 
mK. The k = oT® fittings yield a = 2.90 W/K'*m for Ka 
in 9.95 T and 0.79 W/K"'^m for Kc in 9.75 T, respectively. 
Therefore, the magnon thermal conductivity at Hd can 
be obtained by subtracting the zero-field data from the 
critical-field curves, which gives Km = 2.08T® (W/Km) 
and 0.29T® (W/Km) along the a and c axis, respectively, 
as shown in Figs. 4(c) and 4(d). This nearly perfect T® 
dependence is actually the clearest experimental evidence 
of the magnon ballistic transport in the AF systems till 
now. 

There are several notable details about the field- and 
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temperature-dependencies of k. First, for both field di¬ 
rections, the transition field determined by the peak po¬ 
sitions of k{H) isotherms is nearly temperature inde¬ 
pendent and slightly shifts to lower fields with increas¬ 
ing temperature. This T-dependence of Hd is actu¬ 
ally not the same as the transition field closing the spin 
gap probed by other measurements . Very sim¬ 
ilar phenomenon has been found in another BEC ma¬ 
terial DTNi^ The possible reason is that the magnetic 
excitations may not only transport heat but also scat¬ 
ter phonons. Since these two factors contribute oppo¬ 
sitely to K, the field dependence of k is rather com¬ 
plicated and the maximum could appear at the fields 
slightly different from the phase-transition fields. In 
particular, this discrepancy is stronger at higher tem¬ 
peratures where the scattering between phonons and 
magnons are stronger and the peak becomes broader. 
Second, at the critical fields, k{T) show some devia¬ 
tions of the T-dependence from the exact behav¬ 
ior at both high and low temperatures. At high tem¬ 
perature (> 500 mK), the T-dependence becomes a bit 
weaker, which is a usual phenomenon of heat trans¬ 
port if the scattering between magnetic excitations and 
phonons cannot be neglectedThat is, it is a sim¬ 
ple deviation from the boundary scattering limit. An¬ 
other possible reason for this deviation is a dimensional 
crossover of spin system to the two-dimensional (2D) 
regime, as found, for example, in a 2D antiferromagnet 
RbFe(Mo 04 ) 2 <^ In this regard, there seemed to be no ex¬ 
perimental result indicating such a dimensional crossover 
in IPA-CuClai^i^^^— On the other hand, Ka{T) at 
9.95 T shows stronger T-dependence than T^ at the low¬ 
est temperature regime. So far, the reason is not very 
clear. Probably some decoupling between phonons and 
magnons occurs, like the electron-phonon decoupling in 
high-Tc cupratesj^i^ 

IV. DATA ANALYSIS AND DISCUSSIONS 

Similar to the case of phonons, the ballistic magnon 
thermal conductivity of an isotropic system can be writ¬ 
ten asi“— 

1 2 2 , fkBTY - 2 , 

Km ~ ^ I 7l j ^ 

where Vm is the averaged magnon velocity, and Im is 
the T-independent mean free path of magnons (see Ap¬ 
pendix). It is easy to find that Eq. (jU of the isotropic 
system cannot quantitatively describe the experimental 
results of Km,a and Km,c- Note that at ultra-low temper¬ 
atures, the Im is also determined by the boundary scat¬ 
tering and is actually the same as the mean free path of 
phonons. The calculation using Eq. Q gives the magnon 
velocities of 920 and 1930 m/s along the a and the c axis, 
respectively. They are apparently unreasonable since the 
magnon dispersion is stronger along the a axis (the ladder 
direction) It is simply due to the fact that Eq. (jH) is 


established for the isotropic system and cannot be valid 
for the low-dimensional spin systems. As shown in the 
Appendix, for the anisotropic system the magnon ther¬ 
mal conductivity along a certain direction i is expressed 
as 

1 

^m,i — ^ X 

where Vm,i, lm,i and a' represent the averaged magnon 
velocity along the direction i, the T-independent mean 
free path and the coefficient of frequency density distri¬ 
bution, respectively. Here, a' is written as 

, ^ d6> sin 

/o"^'d0/;/V(0,</)d</ 

while W (9, </) is an assumed velocity-weighting function 
with an exponential form, 

W{9, /)) = e ”0 , (7) 

with vq an adjustable parameter, 9 and (j) the space an¬ 
gles. 

The magnon dispersions along the a and c axis of IPA- 
CuCla have been known from the earlier neutron mea¬ 
surements, which gave Vm,a = 2050 m/s and Vm,c = 790 
m/si^“— The magnon dispersion is known to be very 
weak along the b axis^- Here, we use Eq. m to calcu¬ 
late the Km,a and Km,c with the experimental values of 
Vm,a and Vm,c and adjustable parameters Vm,b and vg- 
According to this formula, these two parameters have 
different impacts on the value of Km,a and Km,c- The 
magnitudes of Km,a and Km,c increase with decreasing 
Vm,b and their ratio (or anisotropy) hardly changes with 
Vm,b- In contrast, Km,a and Km,c decrease with decreasing 
Vq. Furthermore, the ratio of Km,alKm,c increases with 
decreasing vg. Therefore, only with appropriate values 
of Vm,b and vg, one can get calculated k close to the ex¬ 
perimental data for both Km,a and Km,c simultaneously. 
The best calculations of Km,a = 2.27T^ (W/Km) and 
Km,c = 0.26T^ (W/Km), as shown in Figs. 4(c) and 
4(d), are obtained with Vm,b = 150 m/s and vg = 220 
m/s. Note that the parameter Vm,b is about one order 
of magnitude smaller than the velocities along the other 
axes, which is reasonable for IPA-CuCla^ir— Therefore, 
the ballistic magnon heat transport of a low-dimensional 
quantum magnet can be well described by an anisotropic 
formula ([S]). 

Another important result of this work is the peak¬ 
like feature of k{H) at Hd, which becomes clearer and 
sharper with lowering temperature. As already men¬ 
tioned above, this indicates that the spin gap is closed 
only at the critical held and is re-opened in the held- 
induced AF state. Therefore, the lowest excitation in the 
held-induced AF state is non-Goldstone mode. The gap 
size cannot be determined precisely by the present data, 
but we can take a rough estimation from the high-held 


AbTV. 

30’^^ V n ; 


-ks 


,ia 


( 5 ) 
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behavior of k{H). In the high-held state with small gap, 
the magnons can still be easily excited and contribute to 
transporting heat if ksT is not smaller than the gap. As 
shown in Fig. 3, it seems that the n tends to recover its 
zero-held value at high-held limit of || a {H || c ) when 
T < 380 mK (252 mK) but tends to increase at high- 
held limit when T > 520 mK (380 mK). Therefore, the 
gap is estimated to be about 500 mK (« 0.043 meV) and 
300 mK (rs 0.026 meV) for H \\ a and c, respectively. 
Apparently, such small gaps are beyond the resolution 
of the earlier neutron measurements. Similarly, TlCuCla 
was the hrst BEC candidate that showed a gapless Gold- 
stone mode in the held-induced AF state by the neutron 
scattering measurement;^ but was lately found to have 
a small gap of 0.09 meV by the ESR measurement 4^^ 
A characteristic of BEC is the presence of U{1) sym¬ 
metry, which corresponds to the global rotational sym¬ 
metry of the bosonic held phase.— In held-induced 
XY-ordered state, the t/(l) symmetry spontaneously gets 
broken and thus a gapless Goldstone mode is acquired. 
However, the re-opening of the gap is a clear evidence for 
a broken uniaxial symmetry of spin Hamiltonian, which 
rules out a strict description of the magnetic order in 
terms of BEG— Therefore, the heat transport data in¬ 
dicate that the BEC model has limited applicability to 
IPA-CuCls, similar to TlCuCla. The theoretical works 
actually had predicted a general instability of an axi¬ 
ally symmetric magnetic condensate toward a violation of 
this symmetry and the formation of an anisotropy gap at 
H > Hci .— It is related to the presence of anisotropic 
interactions, such as the dipole-dipole coupling, the spin- 
orbital interaction, etc— 4£ 


V. CONCLUSIONS 

In summary, we have studied the ultra-low-T heat 
transport of a spin-gapped compound IPA-CuCla, which 
has been classified to be a candidate of BEC. When the 
gap is closed by the field at Hd, a ballistic magnon 
heat transport is observed. On the other hand, the low-T 
k{H) isotherms show peaks at the critical field, indicating 
that the spin gap is re-opened at H > Hd- Therefore, 
IPA-CuCla seems not to be an ideal BEC prototype sys¬ 
tem. Ultra-low-T thermal conductivity is a very sensitive 
technique to probe the small spin gap, and the validity of 
BEC model for those quantum magnets can be carefully 
examined using this measurement. 
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Appendix: Phenomenological formulas for the heat 
transport of phonons and magnons 

A. The isotropic systems 


According to the Boltzmann equation, the lattice ther¬ 
mal conductivity can be obtained by solving the inte¬ 
gral equation. In an isotropic (or cubic) 3D crystal, the 
phononic thermal conductivity is given bj4r4i42 


K;p - 


1 1 
3 (27r)3 





(Al) 


where Vp^\{q), lp,xiq), Cp^x and fx{q) are the velocity, 
mean free path, specific heat per normal mode and vi¬ 
bration mode distribution function of phonons with po¬ 
larization A (two transverse and one longitudinal) and 
wave vector q, respectively. The specific heat Cp^x is ex¬ 
pressed as 


CpAq) 


dT 


^ujx{q) ■ 


e _ 1 


(A2) 


where ujx{q) is the dispersion for polarization A. 

Usually, the following assumptions are needed to sim¬ 
plify the equation and calculate the result. First, it is 
assumed that phonons of different polarizations have the 
same contribution, which means Vp{q), Ipiq), Cp, f(q) 
and ui{q) are independent of polarization A. Second, 
according to the isotropic Debye approximation, Vp is 
a constant and the dispersion relation is q = uj/vp in 
all directions. Thus, q is independent of direction and 
f{q)dq = dnq'^dq = 47ra;^u“^da;. With these assump¬ 
tions, the thermal conductivity can be simplified as 


1 



r^o/T 


x'^e^ 
(e^ - 1)2 


dx. 


(A3) 


where x = If the temperature is low enough, 

there is only boundary scattering and the mean free path 
is a constant given by the cross-section area, that is. 
Ip = 2^Aj-K. In this case, the integral approaches a 
constant of 47r‘*/I5 and thermal conductivity obeys a T^ 
law, which is formula (jS]) mentioned in the main text. 

This result can also be applied to gapless acoustic 
magnons of a 3D antiferromagnet. Note that when 
the degeneracy of gapped magnon branches is lifted by 
magnetic field, only the lowest branch becomes gapless. 
Hence, the magnon thermal conductivity is 
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I 27r2 

= 3 IVe*'’ 



(A4) 


where Vm a-nd Ijn are the velocity and mean free path of 
magnons. 
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B. The anisotropic systems 

The above formulas for phonons and magnons are valid 
only for the isotropic systems. Usually, the anisotropy 
of phonons is not strong for most 3D crystal lattices and 
therefore formula ([3]) can describe the phonon heat trans¬ 
port in a good approximation. For anisotropic systems, 
these formulas predict that the k is smaller along the di¬ 
rection with stronger dispersion (larger velocity). This 
is apparently unreasonable. For strongly anisotropic sys¬ 
tems, like the magnons in a low-dimensional quantum 
magnet, this problem would be more serious. Here, we 
give some phenomenological results of the Boltzmann 
theory in the anisotropic systems. 

Let us again start with the phonon heat transport. 
Note that with the above assumptions for isotropic 
system fir— phonons only travel along the tempera¬ 
ture gradient direction with a uniform velocity. In the 
isotropic case, we need to consider phonons traveling in 
all directions with component along the temperature gra¬ 
dient. Here, the phonon velocity is a function of space 
angles and a simple function model is an ellipsoid (see 
Fig. 5). Thus, the projections of a random velocity of 
one acoustic branch can be written in parametric form 
as 


Va = a cos a cos 7 , (A5) 

Ub = 6 sin cr cos 7 , (A 6 ) 

'(;c = csin7. (A-7) 

Here, a, b and c are the phonon velocities along the three 
principal crystal axes, a and 7 are parameters of the el¬ 
lipsoid equation and they are represented by space angles 
6 and </> as 


tr = arctan 


itand 
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7 = arctan 


acosa)^ + (fesincr)^ 
ctan(/) 
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Then, the anisotropic sound velocity can be written as a 
function of space angles, 

v{9, <()) = V Va'^{9, (j)) + Vb'^iO, (j)) + Vc’^(9, (j)). (AlO) 

In this anisotropic velocity model, Vp^q and f{q)dq need 
correction. With Eq. (AIO), the averaged phonon veloc¬ 
ity along the a, b and c axis are written as 

d9 v(9,(j)) cos9 sm(j)d(j) 

= - - \/2 J/2 -All 

/;Ad0/;Ad<(> 


Vb = 




(A12) 


(7r/2)/;Ad0/;Ad<^ 


(A13) 


Since v{9, </>) depends on the space angles here, q = 
uj/v{9, 4>) is function of cj, 9 and 4> simultaneously. Thus, 
f{q)dq has a more complex form as 

f{q)dq = Aaoj'^duj, (A14) 


where 


a = 




V ^{9,(j))s\xi(f)d(j). 


(A15) 


Take Eqs. (A2) and (A11-A14) into Eq. (Al), one 
can get the thermal conductivities along direction i for 
phonons and magnons in low temperature limit as 


Kp,^ = ‘^Vp4p,,akB 



and 


1 27r 


i = - X 
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(A16) 


(A17) 


respectively. 

Here, we put forward a further quantitative assump¬ 
tion involving the anisotropy. We consider that phonons 
or magnons prefer to propagate in a direction with 
stronger dispersion and this direction of course con¬ 
tributes more to the average velocity. Thus, we assume a 
velocity-dependent weighting function with an exponen¬ 
tial form, 

•vie,4,) 

W{9,(j)) = e~^, (A18) 


where vq is an adjustable parameter. With this weighting 
function, the averaged phonon velocities in direction a, b 
and c should be written as 


d9 \Y ( 0 ^ (j))v{9, cj)) cos 9 sin (j)dcj) 

f;^^d9f;^^w(9,^)d</> 


(A19) 
















d0 W(ff, (j))v{6, (j)) sin 6 sin (j)d(j) 


(A20) 


Take Eqs. (A2) and (A19-A22) into Eq. (Al), one can 
finally obtain the phonon and magnon thermal conduc¬ 
tivities along direction i (a, b or c) as 


Vr = 


d0 W(0, (/))v(0, (/>) cos 0d(/) 
(7r/2)f;^^d0f;^^W(d,c/>)d</> 


(A21) 


The denominators are normalization coefficients. 

For the same reason, taking the weighting function into 
account, f{q)dq has a similar form of 




and 


^m,i — g ^ h]^ ( ^ ^ ? 


(A24) 


(A25) 


f{q)dq = TT^a'w^dw, (A22) 


and here 


, _ d6) W{0,(l))v ^(0,4>) sin (jydcf) 

r^'d0/;/V(0,</))d<^ 


(A23) 


respectively. 

In the present work, the crystal lattice and also the 
magnetic structure have lower symmetry than the or¬ 
thorhombic one. Using above formulas to fit the experi¬ 
mental data is another approximation. 
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